\(\int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx\) [322]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [B] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [F]
   Giac [A] (verification not implemented)
   Mupad [F(-1)]

Optimal result

Integrand size = 22, antiderivative size = 105 \[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=-\frac {\sqrt {a+c x^2}}{d x}-\frac {\sqrt {c d^2+a e^2} \text {arctanh}\left (\frac {a e-c d x}{\sqrt {c d^2+a e^2} \sqrt {a+c x^2}}\right )}{d^2}+\frac {\sqrt {a} e \text {arctanh}\left (\frac {\sqrt {a+c x^2}}{\sqrt {a}}\right )}{d^2} \]

[Out]

e*arctanh((c*x^2+a)^(1/2)/a^(1/2))*a^(1/2)/d^2-arctanh((-c*d*x+a*e)/(a*e^2+c*d^2)^(1/2)/(c*x^2+a)^(1/2))*(a*e^
2+c*d^2)^(1/2)/d^2-(c*x^2+a)^(1/2)/d/x

Rubi [A] (verified)

Time = 0.11 (sec) , antiderivative size = 105, normalized size of antiderivative = 1.00, number of steps used = 15, number of rules used = 11, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.500, Rules used = {975, 283, 223, 212, 272, 52, 65, 214, 749, 858, 739} \[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=-\frac {\sqrt {a e^2+c d^2} \text {arctanh}\left (\frac {a e-c d x}{\sqrt {a+c x^2} \sqrt {a e^2+c d^2}}\right )}{d^2}+\frac {\sqrt {a} e \text {arctanh}\left (\frac {\sqrt {a+c x^2}}{\sqrt {a}}\right )}{d^2}-\frac {\sqrt {a+c x^2}}{d x} \]

[In]

Int[Sqrt[a + c*x^2]/(x^2*(d + e*x)),x]

[Out]

-(Sqrt[a + c*x^2]/(d*x)) - (Sqrt[c*d^2 + a*e^2]*ArcTanh[(a*e - c*d*x)/(Sqrt[c*d^2 + a*e^2]*Sqrt[a + c*x^2])])/
d^2 + (Sqrt[a]*e*ArcTanh[Sqrt[a + c*x^2]/Sqrt[a]])/d^2

Rule 52

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(a + b*x)^(m + 1)*((c + d*x)^n/(b*(
m + n + 1))), x] + Dist[n*((b*c - a*d)/(b*(m + n + 1))), Int[(a + b*x)^m*(c + d*x)^(n - 1), x], x] /; FreeQ[{a
, b, c, d}, x] && NeQ[b*c - a*d, 0] && GtQ[n, 0] && NeQ[m + n + 1, 0] &&  !(IGtQ[m, 0] && ( !IntegerQ[n] || (G
tQ[m, 0] && LtQ[m - n, 0]))) &&  !ILtQ[m + n + 2, 0] && IntLinearQ[a, b, c, d, m, n, x]

Rule 65

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 223

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Subst[Int[1/(1 - b*x^2), x], x, x/Sqrt[a + b*x^2]] /; FreeQ[{a,
b}, x] &&  !GtQ[a, 0]

Rule 272

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a
+ b*x)^p, x], x, x^n], x] /; FreeQ[{a, b, m, n, p}, x] && IntegerQ[Simplify[(m + 1)/n]]

Rule 283

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[(c*x)^(m + 1)*((a + b*x^n)^p/(c*(m + 1
))), x] - Dist[b*n*(p/(c^n*(m + 1))), Int[(c*x)^(m + n)*(a + b*x^n)^(p - 1), x], x] /; FreeQ[{a, b, c}, x] &&
IGtQ[n, 0] && GtQ[p, 0] && LtQ[m, -1] &&  !ILtQ[(m + n*p + n + 1)/n, 0] && IntBinomialQ[a, b, c, n, m, p, x]

Rule 739

Int[1/(((d_) + (e_.)*(x_))*Sqrt[(a_) + (c_.)*(x_)^2]), x_Symbol] :> -Subst[Int[1/(c*d^2 + a*e^2 - x^2), x], x,
 (a*e - c*d*x)/Sqrt[a + c*x^2]] /; FreeQ[{a, c, d, e}, x]

Rule 749

Int[((d_) + (e_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(d + e*x)^(m + 1)*((a + c*x^2)^p/(e
*(m + 2*p + 1))), x] + Dist[2*(p/(e*(m + 2*p + 1))), Int[(d + e*x)^m*Simp[a*e - c*d*x, x]*(a + c*x^2)^(p - 1),
 x], x] /; FreeQ[{a, c, d, e, m}, x] && NeQ[c*d^2 + a*e^2, 0] && GtQ[p, 0] && NeQ[m + 2*p + 1, 0] && ( !Ration
alQ[m] || LtQ[m, 1]) &&  !ILtQ[m + 2*p, 0] && IntQuadraticQ[a, 0, c, d, e, m, p, x]

Rule 858

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Dist[g/e, Int[(d
+ e*x)^(m + 1)*(a + c*x^2)^p, x], x] + Dist[(e*f - d*g)/e, Int[(d + e*x)^m*(a + c*x^2)^p, x], x] /; FreeQ[{a,
c, d, e, f, g, m, p}, x] && NeQ[c*d^2 + a*e^2, 0] &&  !IGtQ[m, 0]

Rule 975

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))^(n_)*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[ExpandIn
tegrand[(d + e*x)^m*(f + g*x)^n*(a + c*x^2)^p, x], x] /; FreeQ[{a, c, d, e, f, g}, x] && NeQ[e*f - d*g, 0] &&
NeQ[c*d^2 + a*e^2, 0] && (IntegerQ[p] || (ILtQ[m, 0] && ILtQ[n, 0])) &&  !(IGtQ[m, 0] || IGtQ[n, 0])

Rubi steps \begin{align*} \text {integral}& = \int \left (\frac {\sqrt {a+c x^2}}{d x^2}-\frac {e \sqrt {a+c x^2}}{d^2 x}+\frac {e^2 \sqrt {a+c x^2}}{d^2 (d+e x)}\right ) \, dx \\ & = \frac {\int \frac {\sqrt {a+c x^2}}{x^2} \, dx}{d}-\frac {e \int \frac {\sqrt {a+c x^2}}{x} \, dx}{d^2}+\frac {e^2 \int \frac {\sqrt {a+c x^2}}{d+e x} \, dx}{d^2} \\ & = \frac {e \sqrt {a+c x^2}}{d^2}-\frac {\sqrt {a+c x^2}}{d x}+\frac {c \int \frac {1}{\sqrt {a+c x^2}} \, dx}{d}-\frac {e \text {Subst}\left (\int \frac {\sqrt {a+c x}}{x} \, dx,x,x^2\right )}{2 d^2}+\frac {e \int \frac {a e-c d x}{(d+e x) \sqrt {a+c x^2}} \, dx}{d^2} \\ & = -\frac {\sqrt {a+c x^2}}{d x}-\frac {c \int \frac {1}{\sqrt {a+c x^2}} \, dx}{d}+\frac {c \text {Subst}\left (\int \frac {1}{1-c x^2} \, dx,x,\frac {x}{\sqrt {a+c x^2}}\right )}{d}-\frac {(a e) \text {Subst}\left (\int \frac {1}{x \sqrt {a+c x}} \, dx,x,x^2\right )}{2 d^2}+\left (c+\frac {a e^2}{d^2}\right ) \int \frac {1}{(d+e x) \sqrt {a+c x^2}} \, dx \\ & = -\frac {\sqrt {a+c x^2}}{d x}+\frac {\sqrt {c} \tanh ^{-1}\left (\frac {\sqrt {c} x}{\sqrt {a+c x^2}}\right )}{d}-\frac {c \text {Subst}\left (\int \frac {1}{1-c x^2} \, dx,x,\frac {x}{\sqrt {a+c x^2}}\right )}{d}-\frac {(a e) \text {Subst}\left (\int \frac {1}{-\frac {a}{c}+\frac {x^2}{c}} \, dx,x,\sqrt {a+c x^2}\right )}{c d^2}+\left (-c-\frac {a e^2}{d^2}\right ) \text {Subst}\left (\int \frac {1}{c d^2+a e^2-x^2} \, dx,x,\frac {a e-c d x}{\sqrt {a+c x^2}}\right ) \\ & = -\frac {\sqrt {a+c x^2}}{d x}-\frac {\sqrt {c d^2+a e^2} \tanh ^{-1}\left (\frac {a e-c d x}{\sqrt {c d^2+a e^2} \sqrt {a+c x^2}}\right )}{d^2}+\frac {\sqrt {a} e \tanh ^{-1}\left (\frac {\sqrt {a+c x^2}}{\sqrt {a}}\right )}{d^2} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.36 (sec) , antiderivative size = 129, normalized size of antiderivative = 1.23 \[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=\frac {-d \sqrt {a+c x^2}+2 \sqrt {-c d^2-a e^2} x \arctan \left (\frac {\sqrt {-c d^2-a e^2} x}{\sqrt {a} (d+e x)-d \sqrt {a+c x^2}}\right )+\sqrt {a} e x \log (x)-\sqrt {a} e x \log \left (-\sqrt {a}+\sqrt {a+c x^2}\right )}{d^2 x} \]

[In]

Integrate[Sqrt[a + c*x^2]/(x^2*(d + e*x)),x]

[Out]

(-(d*Sqrt[a + c*x^2]) + 2*Sqrt[-(c*d^2) - a*e^2]*x*ArcTan[(Sqrt[-(c*d^2) - a*e^2]*x)/(Sqrt[a]*(d + e*x) - d*Sq
rt[a + c*x^2])] + Sqrt[a]*e*x*Log[x] - Sqrt[a]*e*x*Log[-Sqrt[a] + Sqrt[a + c*x^2]])/(d^2*x)

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(194\) vs. \(2(91)=182\).

Time = 0.38 (sec) , antiderivative size = 195, normalized size of antiderivative = 1.86

method result size
risch \(-\frac {\sqrt {c \,x^{2}+a}}{d x}+\frac {\frac {\sqrt {a}\, e \ln \left (\frac {2 a +2 \sqrt {a}\, \sqrt {c \,x^{2}+a}}{x}\right )}{d}-\frac {\left (e^{2} a +c \,d^{2}\right ) \ln \left (\frac {\frac {2 e^{2} a +2 c \,d^{2}}{e^{2}}-\frac {2 c d \left (x +\frac {d}{e}\right )}{e}+2 \sqrt {\frac {e^{2} a +c \,d^{2}}{e^{2}}}\, \sqrt {\left (x +\frac {d}{e}\right )^{2} c -\frac {2 c d \left (x +\frac {d}{e}\right )}{e}+\frac {e^{2} a +c \,d^{2}}{e^{2}}}}{x +\frac {d}{e}}\right )}{d e \sqrt {\frac {e^{2} a +c \,d^{2}}{e^{2}}}}}{d}\) \(195\)
default \(\frac {-\frac {\left (c \,x^{2}+a \right )^{\frac {3}{2}}}{a x}+\frac {2 c \left (\frac {x \sqrt {c \,x^{2}+a}}{2}+\frac {a \ln \left (x \sqrt {c}+\sqrt {c \,x^{2}+a}\right )}{2 \sqrt {c}}\right )}{a}}{d}-\frac {e \left (\sqrt {c \,x^{2}+a}-\sqrt {a}\, \ln \left (\frac {2 a +2 \sqrt {a}\, \sqrt {c \,x^{2}+a}}{x}\right )\right )}{d^{2}}+\frac {e \left (\sqrt {\left (x +\frac {d}{e}\right )^{2} c -\frac {2 c d \left (x +\frac {d}{e}\right )}{e}+\frac {e^{2} a +c \,d^{2}}{e^{2}}}-\frac {\sqrt {c}\, d \ln \left (\frac {-\frac {c d}{e}+c \left (x +\frac {d}{e}\right )}{\sqrt {c}}+\sqrt {\left (x +\frac {d}{e}\right )^{2} c -\frac {2 c d \left (x +\frac {d}{e}\right )}{e}+\frac {e^{2} a +c \,d^{2}}{e^{2}}}\right )}{e}-\frac {\left (e^{2} a +c \,d^{2}\right ) \ln \left (\frac {\frac {2 e^{2} a +2 c \,d^{2}}{e^{2}}-\frac {2 c d \left (x +\frac {d}{e}\right )}{e}+2 \sqrt {\frac {e^{2} a +c \,d^{2}}{e^{2}}}\, \sqrt {\left (x +\frac {d}{e}\right )^{2} c -\frac {2 c d \left (x +\frac {d}{e}\right )}{e}+\frac {e^{2} a +c \,d^{2}}{e^{2}}}}{x +\frac {d}{e}}\right )}{e^{2} \sqrt {\frac {e^{2} a +c \,d^{2}}{e^{2}}}}\right )}{d^{2}}\) \(370\)

[In]

int((c*x^2+a)^(1/2)/x^2/(e*x+d),x,method=_RETURNVERBOSE)

[Out]

-(c*x^2+a)^(1/2)/d/x+1/d*(1/d*a^(1/2)*e*ln((2*a+2*a^(1/2)*(c*x^2+a)^(1/2))/x)-(a*e^2+c*d^2)/d/e/((a*e^2+c*d^2)
/e^2)^(1/2)*ln((2*(a*e^2+c*d^2)/e^2-2/e*c*d*(x+d/e)+2*((a*e^2+c*d^2)/e^2)^(1/2)*((x+d/e)^2*c-2/e*c*d*(x+d/e)+(
a*e^2+c*d^2)/e^2)^(1/2))/(x+d/e)))

Fricas [A] (verification not implemented)

none

Time = 0.48 (sec) , antiderivative size = 599, normalized size of antiderivative = 5.70 \[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=\left [\frac {\sqrt {a} e x \log \left (-\frac {c x^{2} + 2 \, \sqrt {c x^{2} + a} \sqrt {a} + 2 \, a}{x^{2}}\right ) + \sqrt {c d^{2} + a e^{2}} x \log \left (\frac {2 \, a c d e x - a c d^{2} - 2 \, a^{2} e^{2} - {\left (2 \, c^{2} d^{2} + a c e^{2}\right )} x^{2} - 2 \, \sqrt {c d^{2} + a e^{2}} {\left (c d x - a e\right )} \sqrt {c x^{2} + a}}{e^{2} x^{2} + 2 \, d e x + d^{2}}\right ) - 2 \, \sqrt {c x^{2} + a} d}{2 \, d^{2} x}, \frac {\sqrt {a} e x \log \left (-\frac {c x^{2} + 2 \, \sqrt {c x^{2} + a} \sqrt {a} + 2 \, a}{x^{2}}\right ) - 2 \, \sqrt {-c d^{2} - a e^{2}} x \arctan \left (\frac {\sqrt {-c d^{2} - a e^{2}} {\left (c d x - a e\right )} \sqrt {c x^{2} + a}}{a c d^{2} + a^{2} e^{2} + {\left (c^{2} d^{2} + a c e^{2}\right )} x^{2}}\right ) - 2 \, \sqrt {c x^{2} + a} d}{2 \, d^{2} x}, -\frac {2 \, \sqrt {-a} e x \arctan \left (\frac {\sqrt {-a}}{\sqrt {c x^{2} + a}}\right ) - \sqrt {c d^{2} + a e^{2}} x \log \left (\frac {2 \, a c d e x - a c d^{2} - 2 \, a^{2} e^{2} - {\left (2 \, c^{2} d^{2} + a c e^{2}\right )} x^{2} - 2 \, \sqrt {c d^{2} + a e^{2}} {\left (c d x - a e\right )} \sqrt {c x^{2} + a}}{e^{2} x^{2} + 2 \, d e x + d^{2}}\right ) + 2 \, \sqrt {c x^{2} + a} d}{2 \, d^{2} x}, -\frac {\sqrt {-a} e x \arctan \left (\frac {\sqrt {-a}}{\sqrt {c x^{2} + a}}\right ) + \sqrt {-c d^{2} - a e^{2}} x \arctan \left (\frac {\sqrt {-c d^{2} - a e^{2}} {\left (c d x - a e\right )} \sqrt {c x^{2} + a}}{a c d^{2} + a^{2} e^{2} + {\left (c^{2} d^{2} + a c e^{2}\right )} x^{2}}\right ) + \sqrt {c x^{2} + a} d}{d^{2} x}\right ] \]

[In]

integrate((c*x^2+a)^(1/2)/x^2/(e*x+d),x, algorithm="fricas")

[Out]

[1/2*(sqrt(a)*e*x*log(-(c*x^2 + 2*sqrt(c*x^2 + a)*sqrt(a) + 2*a)/x^2) + sqrt(c*d^2 + a*e^2)*x*log((2*a*c*d*e*x
 - a*c*d^2 - 2*a^2*e^2 - (2*c^2*d^2 + a*c*e^2)*x^2 - 2*sqrt(c*d^2 + a*e^2)*(c*d*x - a*e)*sqrt(c*x^2 + a))/(e^2
*x^2 + 2*d*e*x + d^2)) - 2*sqrt(c*x^2 + a)*d)/(d^2*x), 1/2*(sqrt(a)*e*x*log(-(c*x^2 + 2*sqrt(c*x^2 + a)*sqrt(a
) + 2*a)/x^2) - 2*sqrt(-c*d^2 - a*e^2)*x*arctan(sqrt(-c*d^2 - a*e^2)*(c*d*x - a*e)*sqrt(c*x^2 + a)/(a*c*d^2 +
a^2*e^2 + (c^2*d^2 + a*c*e^2)*x^2)) - 2*sqrt(c*x^2 + a)*d)/(d^2*x), -1/2*(2*sqrt(-a)*e*x*arctan(sqrt(-a)/sqrt(
c*x^2 + a)) - sqrt(c*d^2 + a*e^2)*x*log((2*a*c*d*e*x - a*c*d^2 - 2*a^2*e^2 - (2*c^2*d^2 + a*c*e^2)*x^2 - 2*sqr
t(c*d^2 + a*e^2)*(c*d*x - a*e)*sqrt(c*x^2 + a))/(e^2*x^2 + 2*d*e*x + d^2)) + 2*sqrt(c*x^2 + a)*d)/(d^2*x), -(s
qrt(-a)*e*x*arctan(sqrt(-a)/sqrt(c*x^2 + a)) + sqrt(-c*d^2 - a*e^2)*x*arctan(sqrt(-c*d^2 - a*e^2)*(c*d*x - a*e
)*sqrt(c*x^2 + a)/(a*c*d^2 + a^2*e^2 + (c^2*d^2 + a*c*e^2)*x^2)) + sqrt(c*x^2 + a)*d)/(d^2*x)]

Sympy [F]

\[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=\int \frac {\sqrt {a + c x^{2}}}{x^{2} \left (d + e x\right )}\, dx \]

[In]

integrate((c*x**2+a)**(1/2)/x**2/(e*x+d),x)

[Out]

Integral(sqrt(a + c*x**2)/(x**2*(d + e*x)), x)

Maxima [F]

\[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=\int { \frac {\sqrt {c x^{2} + a}}{{\left (e x + d\right )} x^{2}} \,d x } \]

[In]

integrate((c*x^2+a)^(1/2)/x^2/(e*x+d),x, algorithm="maxima")

[Out]

integrate(sqrt(c*x^2 + a)/((e*x + d)*x^2), x)

Giac [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 146, normalized size of antiderivative = 1.39 \[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=-\frac {2 \, a e \arctan \left (-\frac {\sqrt {c} x - \sqrt {c x^{2} + a}}{\sqrt {-a}}\right )}{\sqrt {-a} d^{2}} + \frac {2 \, a \sqrt {c}}{{\left ({\left (\sqrt {c} x - \sqrt {c x^{2} + a}\right )}^{2} - a\right )} d} + \frac {2 \, {\left (c d^{2} + a e^{2}\right )} \arctan \left (-\frac {{\left (\sqrt {c} x - \sqrt {c x^{2} + a}\right )} e + \sqrt {c} d}{\sqrt {-c d^{2} - a e^{2}}}\right )}{\sqrt {-c d^{2} - a e^{2}} d^{2}} \]

[In]

integrate((c*x^2+a)^(1/2)/x^2/(e*x+d),x, algorithm="giac")

[Out]

-2*a*e*arctan(-(sqrt(c)*x - sqrt(c*x^2 + a))/sqrt(-a))/(sqrt(-a)*d^2) + 2*a*sqrt(c)/(((sqrt(c)*x - sqrt(c*x^2
+ a))^2 - a)*d) + 2*(c*d^2 + a*e^2)*arctan(-((sqrt(c)*x - sqrt(c*x^2 + a))*e + sqrt(c)*d)/sqrt(-c*d^2 - a*e^2)
)/(sqrt(-c*d^2 - a*e^2)*d^2)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {a+c x^2}}{x^2 (d+e x)} \, dx=\int \frac {\sqrt {c\,x^2+a}}{x^2\,\left (d+e\,x\right )} \,d x \]

[In]

int((a + c*x^2)^(1/2)/(x^2*(d + e*x)),x)

[Out]

int((a + c*x^2)^(1/2)/(x^2*(d + e*x)), x)